arXiv:1506.03020vl [cond-mat.quant-gas] 9Jun2015 


Spontaneous formation and non-equilibrium dynamics of a soliton-shaped 

Bose-Einstein condensate in a trap 

Oleg L. Berman, Roman Ya. Kezerashvili, German V. Kolmakov, and Leonid M. Pomirchi 
Physics Department, New York City College of Technology, 

The City University of New York, Brooklyn, NY 11201 
(Dated: June 10, 2015) 

The Bose-stimulated self-organization of a quasi-two dimensional non-equilibrium Bose-Einstein 
condensate in an in-plane potential is proposed. We obtained the solution of the nonlinear, driven- 
dissipative Gross-Pitaevskii equation for a Bose-Einstein condensate trapped in an external asym¬ 
metric parabolic potential within the method of the spectral expansion. We found that, in sharp 
contrast to previous observations, the condensate can spontaneously acquire a soliton-like shape for 
spatially homogenous pumping. This condensate soliton performs oscillatory motion in a parabolic 
trap and, also, can spontaneously rotate. Stability of the condensate soliton in the spatially asym¬ 
metric trap is analyzed. In addition to the nonlinear dynamics of non-equilibrium Bose-Einstein 
condensates of ultra-cold atoms, our findings can be applied to the condensates of quantum well 
excitons and cavity polaritons in semiconductor heterostructure, and to the condensates of photons. 
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I. INTRODUCTION 

Bose-Einstein condensation (BEG) and superfluidity 
are hallmarks of quantum degenerate systems composed 
of interacting bosons. In the BEG state, a substantial 
fraction of particles at low temperatures spontaneously 
occupies the single lowest-energy quantum state thus, 
forming a condensate in the energy space. One of famous 
examples of BEGs is a condensation of ultracold alkali- 
metal atoms at nK range of temperatures The 

advances ultra-low-temperature techniques have already 
lead to the development of various technologies including 
international atomic time keeping, the base of the Global 
Positioning System, which is familiar to everyone. 

Recently, significant progress in solid state physics and 
nanofabrication has enabled to experimentally create a 
class of condensed systems in semiconductor heterostruc¬ 
tures that demonstrate BEG il. In this case, the 
Bose particles are excitons, i.e., electron-hole pairs cou¬ 
pled due to Goulomb attraction in quasi-two-dimensional 
quantum wells (QW), or polaritons, quantum superposi¬ 
tions of excitons and cavity photons (see @ for extensive 
review). The effective mass of these particles is much 
smaller than for their atomic counterpart and it varies 
from a the free electron mass rrie order to ~ 10“^TOe 
depending on the physical realization. As a result, the 
solid-state systems undergo the BEG transition at much 
higher temperatures Tc than the atomic BEG: Tc ranges 
from ^ a few K to ~ 40 K 0,10- Physics of exciton and 
polariton BEGs has already revealed exciting phenomena 
including superfluidity [H, quantized vorticity 0], quan¬ 
tum solitons 0, [l^. a nd a condensed-matter analogue 
of Dirac monopole [ll|. Solid-state BEG is a highly de¬ 
veloping research field due to potential applications in 
quantum and optical computing 112.11311 , nonlinear inter¬ 
ferometry [I^, novel light sources|15j|. and atomtronics 
[T^ . Room temperature BEG has also recently been ob¬ 
served for photons [l3, [il|. In the latter case, the non¬ 


linear interactions between the light quanta, which are 
otherwise linear, are provided by the dye molecules in¬ 
troduced into the microcavity. 

In this paper, we study the dynamics of a trapped 
quasi-two-dimensional BEG in the presence of an exter¬ 
nal source and damping that provides general conditions 
for the Bose-Einstein condensation of particles with fi¬ 
nite lifetime. The reduced dimensionality naturally ap¬ 
pears in the solid-state BEG realizations in planar cavi¬ 
ties 0,10- In trapped atomic BEGs this corresponds to 
the limit case where a characteristic frequency of the trap 
along one direction is much higher to those in two other 
directions M- To capture the experimental conditions 
with spatially asymmetric traps, the condensate dynam¬ 
ics was considered for elliptic traps where the trapping 
potential strengths in two orthogonal directions are not 
the same. In our studies, we perform the simulations 
taking an exciton condensate as a relevant example. 

We report that, in sharp contrast to previous obser¬ 
vations, under certain conditions the condensate sponta¬ 
neously self-localizes in a form of a solitary wave with 
a size smaller than a typical condensate cloud size and 
smaller than the excitation spot size. We found that 
a few types of the soliton-like waves can form, includ¬ 
ing condensate humps and rotating doughnut conden¬ 
sates (rings). Earlier, the solitons propagating on the 
background of a uniform condensate has been observed 
in quasi-two-dimensional solid-state systems 00. In 
those studies the solitons were the condensate perturba¬ 
tions, which have been created artificially by perturbing 
the condensate density by an additional “writing” laser 
beam. In our work, the condensate itself self-organizes 
into a strongly nonuniform soliton-shaped state that be¬ 
haves as a “particle”. This particle travels in a trap 
much like a classical particle that oscillates in an ex¬ 
ternal parabolic potential. Self-organization of the con¬ 
densate into a soliton is caused by the interplay of the 
Bose-stimulated condensate formation and the nonlinear 
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interactions between the Bose particles that is, by the 
universal factors, which are present in any physical real¬ 
ization of a BEC. It was also found that if the eccentric¬ 
ity of the elliptic parabolic trap is high, the solitary wave 
does not form. In the latter case, a conventional fluctuat¬ 
ing condensate, which fills all the energetically accessible 
area in the trap, was formed (20l - l^ . 

The paper is organized as follows. In Sec. [TTl we de¬ 
scribe the model used in the simulations and present the 
method of solution of the nonlinear, driven-dissipative 
Gross-Pitaevskii equation for a BEC trapped in an exter¬ 
nal asymmetric parabolic potential. In Sec. lIIIl we discuss 
our main findings. Our conclusions follow in Sec. IIVI 


II. THEORETICAL FRAMEWORK 

In this section we describe our research methodology 
and present the system parameters, for which the sim¬ 
ulations have been done. To describe the dynamics of 
the dipolar exciton condensate, we utilize the driven- 
dissipative Gross-Pitaevskii equation for the condensate 
wave function. In this work, we take excitons in cou¬ 
pled semiconductor quantum wells, i.e. dipolar exci¬ 
tons, as a physical realization of non-equilibrium Bose- 
Einstein condensate where nonlinearity plays an impor¬ 
tant role d, 0, [23l - l^ . 

The Gross-Pitaevskii equation also captures the dy¬ 
namics of Bose-Einstein condensates of cavity polari- 
tons d] and of ultra-cold atoms [l^. In the latter 
case, the complex terms in the right-hand side of driven- 
dissipative Gross-Pitaevskii equation, Eq. (1) below, cor¬ 
respond to the slow condensate depletion due to the cloud 
evaporation and to the initial injection of relatively hot 
atoms. The effective interaction strength in quasi-two- 
dimensional atomic condensates depends on the conden¬ 
sate density in analogy with the case of dipolar ex¬ 
citon condensates detailed in this Section below. 

The model formulated below captures the realistic de¬ 
tails of pumping in the Bose-Einstein condensed systems. 
Specifically, the particles are injected into the system at 
elevated energies (frequencies) and then relax to the low- 
energy states due to the nonlinear interaction with each 
other. The relevant examples of such systems are men¬ 
tioned in the Introduction. For example, for the dipolar 
exciton condensates, the excitons are created at relatively 
high energies due to coupling of hot electrons and holes 
generated by the external laser radiation [2^. For po- 
laritonic condensates, the scattering into the condensate 
from a thermal bath of non-condensed polaritons occurs 
at the bottleneck energy scale, which significantly exceeds 
the characteristic energies of the particles in the ground 
state Q. In trapped atomic condensates, relatively hot 
atoms with the energies exceeding the ground-state en¬ 
ergy in the trapping potential are initially injected into 
the system. It has already been demonstrates that in 
the case of the dipolar exciton condensates, the account 
for these details results in formation of condensate tur¬ 


bulence under certain conditions [2l|, In this pa¬ 

per we show that this model also predicts formation of 
long-living, soliton-like coherent structures in the con¬ 
densate. In contrast, in the conventional approach (see, 
e.g. Ref. @ for review), the pump rate does not de¬ 
pend on the particle energy thus, the important details 
of the nonequilibrium dynamics of the condensates are 
omitted. 


A. Driven-dissipative Gross-Pitaevskii equation for 
two-dimensional condensates 

At temperatures below the BEC transition tempera¬ 
ture, the dipolar exciton condensate is described by the 
mean-field wave function 'I' = 'I'(r,t), which depends 
on the two-dimensional radius vector in the QW plane 
r = {x, y) and time t. The time evolution of the conden¬ 
sate wave functions in an external trap is captured by 
the Gross-Pitaevskii equation 

ih^ = - —A'l’ + UirW + g'l’l'l’l^ + ihf R- —V- 

ot 2me2, y 2t J 

( 1 ) 

In Eq. o, rrigx is the exciton mass, A is the two- 
dimensional Laplacian operator in the QW plane. The 
parabolic trapping potential for the dipolar excitons is 
U{r) = \ where jx and jy are the poten¬ 

tial strengths in x and y directions, respectively. The last 
term in Eq. o describes creation of the excitons due to 
the interaction with the laser radiation and exciton de¬ 
cay, T is the exciton lifetime. The source term oc R'i >(r , t) 
reflects the fact that the condensate particle creation rate 
~ i9|'I'(r,t)p/9t due to Bose-stimulated scattering into 
the condensate is proportional to the condensate den¬ 
sity |'I'(r,t)p. The effective interaction strength g in 
Eq. ([T|) for the dipolar exciton condensate depends on 
the chemical potential g in the system [2l|, [22[ . In the 
case where the exciton cloud size is much greater than 
the mean exciton separation, the interaction strength is 
g = 27r(e^ZI^/i/e^)^'^^, where e is the electron charge, D 
is the interwell distance, e is the dielectric constant of the 
material in the gap between two quantum wells [2l|, . 

To study the exciton condensate dynamics, we numeri¬ 
cally integrate Eq. (HD. There are numbers of approaches 
for the numerical integration of Gross-Pitaevskii-type 
equations. One of the approaches is in solving Eq. (HD 
in r-space by discretizing it using a Grank-Nicholson fi¬ 
nite difference scheme j29l - l3ll |. In Refs. IH and [13 the 
ground-state wave function of a trapped BEC was found 
by the direct minimization of the energy functional for 
the Gross-Pitaevskii equation. The convergence of vari¬ 
ous methods has been studied in Refs. [13 and The 
stability and time-evolution of solutions of the driven- 
dissipative multidimensional Gross-Pitaevskii equation 
has recently been analyzed numerically in Ref. 1361 The 
spectral representation of the Gross-Pitaevskii equation 
has been considered for atomic BEC condensates with- 



3 


out a trap [s^- In that case, the condensate wave func¬ 
tion was expanded using the plane-wave basis. Compre¬ 
hensive reviews for these methods are given in Refs. [H 
and [ 3 ^ 

In our approach, we use the spectral representation for 
the condensate wave function '!'(/•, t) by expanding it in 
terms of basis functions of the exactly solvable station¬ 
ary eigenvalue problem and for the time-dependent coef¬ 
ficients of the expansion obtain the system of the first- 
order differential equations that we solve numerically. 

To utilize our approach let us consider the linear Her- 
mitian part of the Hamiltonian of Eq. © 


After the substitution of the expansion 0 into Eq. 
© one obtains the following system of the first-order 
differential equations for the coefficients An{t), 

= huJnAn{t) + ihy^RnmAmjt) - i^An{t) 
m 

+ g 'y ] WnmpqAmit)Ap(t)A^(t), (8) 

m,p,q 

where 

Wnmpq = ^QWnxm.xPxqx'^nymyPyqy j 


2mex 



( 2 ) 


which is the Hamiltonian of a two-dimensional asymmet¬ 
ric harmonic oscillator. This Hamiltonian enables the 
variable separation and factorization of the wave func¬ 
tions [ 4 ^ . 


Ho'fpnx{x)'ipny{y) = Er,y,ny'llJnx{x)'llJuy{y), (3) 


where + E^l, E^i = and E^l = 

fiLOnl determine the oscillator eigen frequencies and 
LOnl y—direction, respectively. The functions 

and '>pnl{y) are the eigenfunctions of a classical 
one-dimensional harmonic oscillator problem that obey 
the time-independent Schrddinger equation 

- (4) 

where ^ labels the x and y variables. The solution of 
Eq. © is [13 


i’tx (a;) 
i^uiiy) 


1 1 

sj ria;! 
1 1 




(5) 

,( 6 ) 


where -b 1/2), = 

-b 1/2), Ux = 0,1, 2,..., Uy = 0,1, 2,..., 

Hn are the Hermite polynomials, and Xx and Xy charac¬ 
terize the oscillator length scales in the x and y directions, 
respectively. 

To capture the presence of the trap now let us ex¬ 
pand the condensate wave function 4'(r,t) in terms of 
the eigenfunctions ipnir) = 


='^Ar,{t)-ipnir), (7) 

n 


/ oo 

dUniOfmiOfpiOfqiO, 

-OO 

and 

vV/22Ihi! 

The matrix elements of the source operator R in Eq. 
(fTTl) are 


Rnm — / dv ll)xi{v')R’ll^rn{x')' 


(9) 


We consider the case where the excitons are created by an 
external homogenous source in a given range of energies 
huj 2 )j where wi and uj 2 are the boundary frequencies 
of the excitation frequency range. Thus, we assume the 
following form for the matrix elements ©, 

Enm — (19) 


where A(a;n) = i?o if A A <-^2 and A{uin) = 0 
otherwise. The Ro constant characterizes the intensity 
of the exciton source. In most simulations, it was set 
equal Rq = O.ISwq where wq = ( 70 /wea;)^^^ is the oscil¬ 
latory unit of frequency. Additionally, the simulations 
have been performed for Rq = 0.25a;o to demonstrate the 
transition to turbulence at elevated pump rates. In the 
simulations, we set wi = 5wo and 0^2 = Twq. 

We consider Eqs. © and © in the interaction repre¬ 
sentation by separating the time dependence for the lin¬ 
earized equation, An{t) = an(t)e“®‘^"‘, where ujn is the 
eigenfrequency of the mode n. In this representation, 
Eq. © reads 

\Rnn 2 ) 

+ g'^Wnmpq amit)ap(t)a*(t)e^^‘^\ (II) 

m.p.q 

where Wnmpq nre the matrix elements of the dipolar exci¬ 
ton interaction. Aw = Wq-bw^ —Wm —Wp is the frequency 
detuning, and a star stands for the complex conjugate. 

The characteristic length scales of the linearized prob¬ 
lem are 


A 


X 



m'jy 


where An{t) are the time-dependent coefficients of the 
expansion, and n = (hx, Uy) is a two-dimensional integer 
index. 


1/4 


( 12 ) 
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We assume that the trap is elongated in x direction thus, 
Aa: > The characteristic trapping potential length 
scale is Ap = y/XxXy and the mean trapping potential 
strength is 70 = y/'jx'jy- To characterize the spatial 
asymmetry of the trapping potential U{r) we introduce 
its eccentricity as follows, 


e = 


XI-XI 


(13) 


While the dependence of the condensate dynamics is 
studied for different eccentricities e, the total area ac¬ 
cessible for the exciton cloud S = ttAq in the trap U{r) 
is fixed to keep the average condensate density constant 
at a given particle creation rate. 

The trap parameters 'jx and 7 y, which represent the 
confinement strength, are expressed through the eccen¬ 
tricity (fT^ as 


72 . =7o(l-e^), = (14) 


The characteristic length scales in the x and y directions 
are 


Ax — 


Ao 


(1- £2)1/4’ 


Ay = Ao(l-£2)1/4. (15) 


The eigenfrequency of the oscillatory mode n = {nx,ny) 
in a trap with the eccentricity e is 


^Tl 


70 

'^ex 



- e^(nx + 1 / 2 ) 



(riy -t- 1/2) 


(16) 


The initial conditions at t = 0 were set as a Rayleigh- 
Jeans-like thermal distribution 


an(0) 


ksT 

fj,0 + tux)n 



(17) 


with the chemical potential = ^p, random phases 
(j)n and the temperature T = O.lhujo/kB, where ks is the 
Boltzmann constant. 


B. Simulation parameters and the integration 
method 

As we stated above, we consider dipolar excitons 
with spatially separated electrons and holes as the main 
example of non-equilibrium Bose-Einstein condensation 
[ 1 , 0 , [2a - [^ . Taking GaAs heterostructures as a relevant 
example of such system, we set the dielectric constant 
equal e = 13. The simulations have been done for yp = 50 
eV/cm^, the exciton mass mex = 0.22me where nie is 
the free electron mass, the interwell distance D = 4.2 
nm, and the exciton lifetime r = 100 ns, which are the 


representative parameters for the experiments with the 
dipolar excitons in GaAs coupled QWs [13,113, H^]. In 
the simulations, we express the spatial coordinates and 
time in the oscillatory units of length Ax = Ap = 0.9 
/im and of time At = (jnexl^oY^^ — 1-6 ns. The unit of 
frequency in the simulations is wp = At~^ = 6.3 x 10® 
S“4. 

As it is detailed above, we consider the dipolar exciton 
condensates under the conditions of the experiments in 
Refs, i, 0 mu where a rarefied exciton gas with the 
density of n 2 D ^ 10® — lO^® cm“^ is formed in coupled 
QWs with the interwell separation D > O.Sqb, where 
gb = eh^/2mre'^ Ri 7 nm is the two-dimensional exci¬ 
ton Bohr radius and is the exciton reduced mass. 
Such dilute electron-hole system {aBn^jj ^ 4) can 
be described as a weakly-non-ideal Bose-gas of excitons 
[U, m|. As it was shown in Refs, m and m, in this 
regime a gas of interacting dipolar excitons can be consid¬ 
ered as a quantum fluid and the probability of biexciton 
formation is negligibly small. 

We numerically integrated Eq. (EH) with the 4th or¬ 
der Runge-Kutta scheme with the time step of lO^^At. 
To obtain the relative accuracy better than 10“® for the 
condensate wave function and the relative accuracy of 
4 X 10“® for the total number of excitons in the con¬ 
densate, we use Nx x Ny = 256 basis functions in the 
expansion © where 0 < < A^x, 0 < riy < Ny, and 

Nx=Ny = 16 HI. 

In the present work, we utilized Graphics Processing 
Unit (GPU) NVIDIA Tesla K20m to solve numeri¬ 
cally Eqs. EH- The use of GPUs enabled us to perform 
the simulations within a practically reasonable amount 
of time and to achieve the required accuracy. In the 
GPU implementation, we separated the real and imagi¬ 
nary parts in Eqs. EH) and thus, obtained a set of real 
equations. The interaction amplitudes Wnmpq and the 
coefficients a„(0) were initialized on the host (proces¬ 
sor) and then, were transferred from the host memory 
to the device (GPU) memory prior to the start of the 
integration. In the GPU implementation of the code, we 
used a two-dimensional (2D) computational grid consist¬ 
ing of 2D blocks, which in their turn consisted of a set of 
threads that worked in parallel. We divided the problem 
size, Nx X Ny, in such a way that each thread on GPU 
operated on one element of the resultant grid that is, up¬ 
dated the variables at given {nx,ny) at the current time 
step. 

It follows from Eqs. EH that for each n, the nonlin¬ 
ear interaction term involves a summation over six one¬ 
dimensional indices. Given the 2D computational grid, if 
each thread in a block would access the GPU memory to 
read the data, the performance of the kernel would suffer 
from the memory bandwidth bottleneck in addition to 
the read latency. To effectively utilize the memory band¬ 
width and avoid redundant read accesses, we performed 
the following optimizations. The GPU kernel has been 
divided into three sub-kernels, which performed partial 
summations in the nonlinear term of Eqs. (HID . The 
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first sub-kernel produced the summation over the two- 
dimensional index q. The output of this kernel has been 
stored in the global memory on the device and then, has 
been used by the second partial summation sub-kernel 
to produce the partial sum over the index p. The latter 
partial sum was used by the third sub-kernel to produce 
the hnal sum in the nonlinear term in Eqs. m- It was 
found that splitting the summation kernel into three sub¬ 
kernels results in lower execution time than performing 
the complete summation in a single kernel. 

We used a serial code for the same parameter set to 
benchmark the GPU code. In the serial code, the explicit 
complex representation of Eqs. (HU has been utilized. 
We found that the simulation for a single run required 
~ 88 hours on a GPU and showed ^ lOx acceleration 
compared to the serial version of the code. 

III. RESULTS AND DISCUSSION 



FIG. 1. Dependence of the total number of dipolar excitons 
in the condensate on time in a trap with 70 = 50 eV/cm^ for 
the trap eccentricity e = 0.2. The numerical unit of time At 
is specified in the text. 


In this section we present the results of our studies. 
Eirst, we consider spontaneous formation of the dipo¬ 
lar exciton condensate patterns and their dynamics in a 
trapping potential. Then, we analyze the details of the 
pattern rotation. Finally, we consider the effect of the 
asymmetry of the trapping potential and compare our 
findings with the existing results. 

A. Formation of a soliton-like condensate 


pattern is formed. Below we refer to such structures as 
condensate solitary waves, or condensate solitons. 

It is also seen in Fig.[5J the pattern travels in the trap¬ 
ping potential. To characterize the condensate pattern 
motion, we show in Fig. [3K,b the dependence on time of 
the X and y coordinates of the center of mass of the dipo¬ 
lar exciton condensate for the data presented in Fig.[2^-c. 
The center-of-mass coordinate remit) = {xcm{t),ycm(t)) 
was calculated from the condensate wave function as fol¬ 
lows 


To study the dynamics of the dipolar exciton conden¬ 
sate, we determined the condensate density distribution 
at different moments of time by integrating Eq. (HD with 
the source term HHD and initial conditions HZD- In our 
studies, we performed six independent runs for the trap 
eccentricity e = 0.2 and three independent runs for all 
other values of the eccentricity in the range 0 < e < 0.8. 
Fig. [T] demonstrates the dependence of the total num¬ 
ber of dipolar excitons in the condensate as a function of 
time, 


iV(t)=y'dV|d/(r,t)|2, (18) 

obtained in the simulations for the trap eccentricity e = 
0.2. It is seen that initially the number of excitons grows 
with time and then, saturates aX t > ts = 70At Ri 1.1 x 
10“^ s. At t ^ ts the system comes to a steady state 
where the total number of dipolar excitons only slowly 
varies with time. 

To characterize the dipolar exciton condensate dynam¬ 
ics in the steady state at t > G, we studied the time evo¬ 
lution of the condensate density in the trap. Fig.[2]shows 
the graphical output from the simulations made for the 
eccentricity e = 0.2 (upper row) and e = 0 (lower row). 
The results in Fig.[2^-c are shown for the same data with 
e = 0.2 as in Fig.[T]above. It is seen in Fig.[2]that in both 
cases e yf 0 and e = 0, a spatially-localized condensate 


„ /d^r'|^'(r,t)|V 


(19) 


It is seen from Fig.[3^,b that the position of the center of 
mass of the pattern oscillates with time. At the moment 
t ~ ISOAf, the oscillations reach a steady state, in which 
their amplitude and period only slowly varies with time. 
As follows from Fig. [3^,b, the period of the oscillation of 
the pattern in the trap is ^ 6.3At Ri 1.0 x 10“® s. The 
trajectory of the pattern can also be viewed as a Lissajous 
figure for the condensate center of mass in a parametric 
plot in the {xcmit),ycm{t)) plane, as it is presented in 
Fig. [3t. The distance between two farthest points of the 
trajectory for the condensate center of mass estimated 
from Fig. [5^-c is ^ 4.8Aa; r; 4.3 ^m. 

To characterize the internal structure of the soliton- 
like pattern, we calcnlated the time- and angle-averaged 
radial density distribution function for the condensate in 
the center-of-mass frame. 


D{p) = difipl'i’iremit) + p,t)\^'^ . 

( 20 ) 

In Eq. d^HD, (G — ti) is the time-averaging interval, the 
integral in the round brackets is taken over the direction 
ipp of the relative radius-vector p = r — remit), and we 
denoted p = \p\. The time-averaging interval was bound 
by the moments ti = 140At and t 2 = 200At. The results 
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FIG. 2. (Color online) Formation of a traveling soliton-like dipolar-exciton pattern in a trap with the eccentricity e — 0.2 
(upper row) and e = 0 (lower row). The hgures show the exciton condensate density n{r,t) = |^'(r,t)p in three subsequent 
moments. Upper row: (a) t = 139.5At, (b) t = 142.5At, and (c) t = 145.5At after the external drive was turned on. Lower 
row: (d) t = 143.OAt, (e) t = 146.5At, and (f) t = 149.5At. The vertical bars show the exciton condensate density. The 
numerical unit of length Ax is specified in the text. 


of the calculations for three independent runs for e = 
0 and £ = 0.2 are shown Fig. SI It is seen that the 
characteristic radius po of the pattern defined as a width 
at the half height for D{p) varies from ^ 2Aa: to 2.5Aa:. 
This conclusion was validated for all independent runs 
performed under the same conditions. The comparison of 
Figs.[3^-c and Fig. SI shows that the traveled distance for 
the soliton-like pattern during its oscillation in the trap 
is of the order of or greater than twice the radius of the 
pattern. In other words, a soliton-like pattern performs 
large-amplitude oscillatory motion in a parabolic trap. 

It is also seen in Fig. SI that two types of patterns are 
formed. In the pattern of the first kind, the maximum of 
the exciton condensate density is positioned at its center 
'T'cm (see open circles in Fig. SI). This corresponds to a 
hump-like soliton structure. In the pattern of the second 
kind, the condensate density has a minimum at the cen¬ 
ter of the pattern (open triangles and filled triangles in 
Fig. SI and therefore, it corresponds to a ring-like struc¬ 
ture. We also found that the structures of both kinds can 
be formed for the same eccentricity (cf. open circles and 
open triangles for £ = 0 in Fig.|4|). It is worth noting that 
the runs are independent since each run was set based on 
random phases in the initial conditions for the conden¬ 
sate wave function at f = 0 as described in Sec. Ill Al The 
motion of a ring structure is shown in the first row of 
Fig. [21 for which the minimum of the condensate den¬ 
sity is seen as a bright spot at the center of the pattern, 
whereas a hump-like structure is shown in the second row 
of the same figure. 

We infer that the solitary-wave condensate formation 


is caused by the interplay of the two factors: (i) Bose- 
stimulated scattering of excitons into the condensate and 
(ii) the exciton-exciton repulsion in the condensate. The 
first factor results in the increased probability oc Ui of the 
condensate growth in the state i in the region already oc¬ 
cupied by Tii particles [l^ . This effect is captured by the 
source term in Eq. m as explained in Sec. Ill Al It leads 
to the exponential-like increase of the particle occupation 
numbers in the occupied states. However, if the spatial 
density of the condensate is increased, the mutual re¬ 
pulsion of the condensate particles tends to decrease the 
density and results in “spreading” of the particles over 
the system. The nonlinear term in Eq. o is responsible 
for this effect. In the presence of the continuous source of 
the particles (i.e., external laser radiation), this is the fi¬ 
nite lifetime of dipolar excitons that limits the growth of 
the total number of the particles in the condensate. The 
approach where the particles are injected into the system 
at a given energy scale has earlier been used to describe 
the dynamics of the Bose-Einstein condensation of ultra¬ 
cold gases [U m, (see Ref. [S^ for review). Specif¬ 
ically, in that scenario, the injected hot atoms relaxed 
due to multiple collisions toward the condensate state. 
The condensate depletion due to atom evaporation was 
captured by introducing the effective condensate decay 
rate. In contrast, for the case of dipolar exciton conden¬ 
sate considered in this paper, the depletion rate is set 
the same for all states to capture the fact that the dipo¬ 
lar exciton lifetime only weakly depends on the exciton 
energy. 
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FIG. 3. Motion of (a) x- and (b) ^-coordinates of the center 
of mass of the condensate pattern. Figure (c) shows motion 
of the center of mass of the condensate in the {x, y) plane for 
t > ISOAf, where the pattern exhibits stable periodic oscilla¬ 
tions in the parabolic trap. The eccentricity of the trapping 
potential is e = 0.2. The curves in figures a-c are shown for 
the data in Fig. [^-c. 


B. Ring condensate rotation 

We assume that formation of the ring structures is 
caused by spontaneous rotation of the non-equilibrium 
condensate. It is known that for the Bose-condensed, 
superfluid systems two types of rotation are possible: 
(i) vortex rotation and (ii) solid-body-like rotation [^ . 
The vortex rotation is realized as quantized vortices in 
the condensate. The latter have a velocity singularity 
and zero condensate density at the vortex core |5l| . In 
addition to classical experiments with superfluid He-II 
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FIG. 4. Time-averaged radial density distribution dSOl) of the 
dipolar exciton condensate soliton as a function of the radial 
distance p in the center-of-mass frame for three independent 
runs with the eccentricity e = 0 (open circles and open trian¬ 
gles) and for e = 0.2 (filled triangles). The averaging is made 
for the time period between the moments ti = lOOAt and 
t 2 = 200At. Data for e = 0 (open circles) and e = 0.2 (filled 
triangles) are calculated for the data shown in Fig. [2] in the 
upper and lower rows, respectively. Points show the result of 
the calculation, curves are shown to guide the eye. 


(see, e.g., a review in Ref. Is^) . quantized vortices have 
recently been observed in Bose-Einstein condensates of 
ultra-cold atom s 15311 and in condensates of exciton po- 
laritons [13, l55j|. For the solid-body rotation, the 
condensate density does not mandatory turn to zero at 
the rotation center [l^. Below we demonstrate that the 
solid-body rotation is realized for the soliton-like patterns 
described above. 

To analyze the rotation of the pattern we, first, calcu¬ 
late the angular momentum of the dipolar exciton con¬ 
densate as a function of time 

M{t) = [ (21) 


where 


M = [p X p] 

is the angular momentum operator for the condensate 
in the center-of-mass frame, p = —ihVp = is 

the linear momentum operator, and 2 ; marks the direc¬ 
tion perpendicular to the QW plane. Fig. [5^ demon¬ 
strates the time dependence of the angular momentum 
of the ring pattern shown in Fig. [2^-c. We see in Fig. [SK 
that the angular momentum of the ring pattern fluctu¬ 
ates with time. We note that the angular momentum 
fluctuations are natural for the dipolar exciton conden¬ 
sate since the system is driven by an external pump: ac¬ 
cording to Eq. dD) both the source and damping of the 
excitons are present. Additionally, the ring condensate is 
formed in the parabolic trap that breaks the rotational 
symmetry in the moving, off-center frame of reference. 
Therefore, the angular momentum of the system is not 
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FIG. 5. (a) Dependence of the angular momentum of the ring 
structure in Fig. [2^-c on time. Points show the results of the 
calculation, the line segments between the points are shown 
to guide the eye. The momentum is expressed in the units 
of h. (b) The time-dependence of the probabilities \cm{t)\^ 
of having the momentum m. (c) Time-averaged probabilities, 
(|cm(t)|^), for various momenta m. The graphs in a-c are 
plotted for the data shown in Fig. [2K-c 


conserved. However, it is seen in Fig. [S^ that at the mo¬ 
ment of time t ~ 140At the angular momentum relaxes 
to a steady state and it further fluctuates about the mean 
value M Ri 0.8. 

To more fully characterize the pattern rotation, we cal¬ 
culated the time dependence of the probabilities |cm(t)P 
of having the specific angular momentum m. Here, Cm{t) 
are the normalized components the expansion of the ex- 
citon condensate wave function over the basis functions 
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FIG. 6. (Color online) Exciton density distribution in a trap 
with the eccentricity e = 0.8 for the moments (a) t = 139.5At, 
(b) t = 142.5At, and (c) t = 145.5At. The vertical bars show 
the exciton condensate density n{x,y). It is seen that a trav¬ 
eling exciton pattern does not emerge for a large eccentricity 
£, and a turbulent condensate [U, is formed instead. 


with a fixed momentum m in the center-of-mass frame, 

/ d(ppdpp^>{r,t)^*^{ipp) 

Cm{t) = — 

(J d(ppdpp\'ii{r,t)\A 

where in Eq. (1^ one has r = r^m + P, I'cm is defined in 
Eq. ([19]), is the eigenfunction of a 

system with a given angular momentum m. The compo¬ 
nents (Ea are normalized so that M{t) = J2m 
Fig-Et* shows the results of the calculations for |cm(t)P 
for the pattern in Fig. Specifically, the components 
with the angular momenta m = 1, 5, and 9 are plotted. 
It is seen that the probabilities |cm(t)P strongly fluctuate 
with time. These fluctuations lead to the fluctuations of 
the angular momentum M{t) of the pattern seen in Fig. 
ISi. The fluctuations are caused by the energy redistribu¬ 
tion between the components with different m due to the 
nonlinear interactions in the condensate as well as due to 
the external pump. 

Figure jSJ: shows the time-averaged probabilities 

(|c^P) = -^ (23) 

*2 — Jti 
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for the components with —20 < m < 20. The prob¬ 
abilities in Eq. (1231) are averaged over the time inter¬ 
val between ti = 140At and t 2 = 200At. It is seen in 
Fig. [5j: that despite the averaged angular momentum of 
the system oscillates around the mean value M ^ 0.8, 
the rotational motion of the pattern is more complex 
and includes the components with the momenta reach¬ 
ing \m\ ^ 10. Moreover, both positive and negative m 
contribute to the total angular momentum therefore, ro¬ 
tations in both directions are present in the ring pattern 
at the same time. If a vortex-like rotation were realized 
in the system, the condensate density inevitably turned 
to zero at the center of the vortex, as is described above. 
Additionally, the vortices with |m| > 1 are dynamically 
unstable and tend to decay to elementary vortices with 
m = 1 [i^. However, in our case, the condensate den¬ 
sity does not turn to zero at any point inside the pattern 
for both hill-like and ring-like structures (see Fig. a and 
its spontaneous rotation includes the components with 
|m| > 1 as is demonstrated in Fig. [SJ:. Therefore, the 
solid-body-like rotation is realized in the case of the dipo¬ 
lar exciton condensate patterns. We emphasize that this 
average solid-body-like rotation is a superposition of a 
few interacting rotating flows of the condensate in the 
pattern, which interference results in the fluctuations of 
the condensate density of the pattern, as is clearly seen 
in Fig. [2] above. 




C. Effects of the eccentricity and pumping rate 

We also studied the effects of the eccentricity e of the 
trapping potential on the condensate dynamics. It was 
found that for symmetric or slightly asymmetric traps 
with the eccentricity £ < 0.5 the oscillating soliton pat¬ 
terns were formed similarly to those described in the pre¬ 
vious sections. However, for relatively large eccentricities 
£ > 0.5 formation of the oscillating patterns was not ob¬ 
served. The result of the simulations for £ = 0.8 is shown 
in Fig. [6] The absence of the oscillating patterns for large 
£ can be understood as follows. For a large eccentricity, 
the width of the spatial domain in the y-direction en¬ 
ergetically accessible for the condensate in the trapping 
potential becomes comparable with or smaller than the 
radius of the pattern po and the pattern is squeezed in 
the 2 /—direction by the trap. Due to that, the pattern 
becomes unstable and the exciton cloud filling all ener¬ 
getically accessible spatial domain in the trap is formed, 
as is seen in Fig.[6l In this case, the observed condensate 
dynamics are similar to a turbulent exciton condensate 
described in Refs. [2l| and [22|, where nonlinear fluctua¬ 
tions of the cloud density are observed rather than the 
formation of a traveling, spatially localized soliton-like 
condensate structure. 

It was also found that the soliton-like condensate is 
not formed even for e <C 1 if the particle creation rate 
i?o is higher than a certain value. In the latter case, a 
turbulent condensate was observed in the simulations, in 



FIG. 7. (Color online) Evolution of the condensate dy¬ 
namics from a traveling soliton (a,b) to a turbulent conden¬ 
sate (c,d) with the increase of the particle creation rate from 
Ro = 0.15ti;o to = 0.25a;o. The eccentricity of the trap is 
£ = 0, time t is shown on the plot. 

agreement with the previous works Refs. and [22|. In 
this turbulent state, a relatively large fluctuating cloud 
is formed at the center of the trap. The transition from a 
soliton-shaped to turbulent condensate with the increase 
of the particle creation rate Rq is shown in Fig. [71 


D. Comparison of our findings with the existing 
results 

The main object to which our findings can be applied is 
a non-equilibrium solid-state condensate, e.g. a conden¬ 
sate of excitons, bound electron-hole pairs, in semicon- 
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ductor heterostructure. Since the lifetime of excitons in 
a single semiconductor QW is relatively short (~ a few 
ns) [56 


the experiments with exciton condensates are 
mostly focused on so-called dipolar (or, indirect) exci¬ 
tons in coupled QWs, where the positively charged holes 
and negatively charged electrons are located in different 
QWs, which are separated by a nanometer-wide semicon¬ 
ductor or dielectric barrier Bin [ 53 . The dipolar ex¬ 
citons formed by spatially separated charges in coupled 
QWs are characterized by the relatively long exciton life¬ 
time compared to excitons in a single QW due to small 
electron-hole recombination rates suppressed by the bar¬ 
rier between QWs. It is worth noting that various ex¬ 
periments have revealed rich collective dynamics of dipo¬ 
lar excitons in coupled QWs and have demonstrated the 
existence of different phases in an electron-hole bilayer 
system. The experimental works on excitonic phases in 
coupled QWs were reviewed in Ref.[5^ The experimental 
progress toward probing the ground state of an electron- 
hole bilayer by low-temperature transport was reviewed 
in Ref. [^. The recent progress in the theoretical and 
experimental developments in the studies of a dipolar 
exciton condensate in coupled QWs was reviewed in Ref. 


Owing to nonzero average electric dipole moment of 
dipolar excitons, their interactions are long-ranged and 
decrease with the distance as oc r~^. The presence of 
the long-range interactions can signihcantly modify the 
condensate dynamics under certain regimes M- In par¬ 
ticular, in spatially homogenous two-dimensional systems 
the logarithmic corrections to the chemical potential of 
the system become essential [^. However, it was demon¬ 
strated [ 21 I, [23 that under usual experiment conditions 
for the dipolar exciton condensation in coupled QWs, 
the most significant contribution to the exciton-exciton 
dynamics arises from short-range scattering. In effect, 
the dynamics of the dipolar exciton condensate can effec¬ 
tively be described as the local one, where, however, the 
effective interaction strength for the condensate particles 
pairwise interactions is a function of the chemical poten¬ 
tial in the system. This quasi-local model has recently 
been applied to the studied of the nonlinear evolution of 
dipolar exciton condensates in a radially-symmetric trap 

11 [23. 

The observed self-organization of the dipolar exciton 
condensates may be similar to the dissipative soliton for¬ 
mation known for open nonlinear optical systems, lasers, 
and chemical systems [^. Solitons dynamics described 
by one-dimensional conservative Gross-Pitaevskii equa¬ 
tion without the source and decay terms (also known as 
nonlinear Schrodinger equation) was studied in details for 
the case where the trapping potential is absent (see Refs, 
m and m for review). However, in our case, the sys¬ 
tem is two-dimensional that makes inapplicable the gen¬ 
eral analytical approaches developed for the integrable 
one-dimensional nonlinear systems . An extended di¬ 
mensionality of the systems compared to Refs. IH andlH 
results in the presence of additional degrees of freedom, 


in particular, in the possibility for the soliton to rotate. 

It is also worth noting that the soliton pro pag ation 
has recently been observed in Refs. [3 [l3, and IbSl for a 
condensate of polaritons, which are a quantum super¬ 
position of the QW excitons and cavity photons Q- In 
Ref. [3 the solitons propagate as a perturbation of the 
steady-state condensate density. In Refs. [13 and [03 the 
solitons were created by a tightly focused “writing” laser 
beam and then, freely expanded or traveled in the mi¬ 
crocavity. In our case, the exciton condensate spon¬ 
taneously forms a traveling soliton under spatially ho¬ 
mogenous pumping that is, inside the pumping spot. 
Formation of inhomogeneous polariton condensate pat¬ 
terns have also been considered in Refs, [gl and [63 In 
those works, the inhomogeneities in the condensate den¬ 
sity were caused by the interactions of polaritons with 
a cloud of non-condensed excitons and with structural 
defects in the cavity. The ring pattern formation was 
earlier demonstrated for the dipolar excitons in coupled 
QWs and multiple QWs in previous experimental and 
theoretical works [3, [SJ, [53 [13 [13 ■ In spatially resolved 
measurements, the photoluminescence of the dipolar ex¬ 
citon pattern was seen as a central laser excitation spot 
surrounded by a sharp ring-shaped region with a diame¬ 
ter much larger than that for the spot and with large dark 
regions separating the two. In the coupled QWs experi¬ 
ments the voltage was applied perpendicular to the wells, 
causing the bands to tilt and the electrons and holes to 
separate. Both the lifetime and the energy of dipolar 
excitons with spatially separated electrons and holes in 
different QWs can be tuned by the electric field. The life¬ 
time of quasiparticles in those systems was increased up 
to microseconds due to the charge separation in different 
QWs [ 73 . In the cases mentioned above the formation 
of the external rings was caused by the recombination 
of in-plane spatially separated electrons and holes. The 
electrically injected carriers created by the external gate 
voltage recombine with optically injected carriers of the 
opposite sign [H, [^ [73 ■ In the framework of the micro¬ 
scopic approach formulated in terms of coupled nonlinear 
equations for the diffusion, thermalization and optical de¬ 
cay of the particles, the formation of the inner ring was 
explained by the fact that in the optically-pumped area 
the exciton temperature is much larger than the lattice 
temperature. As a result, the recombination decay of 
excitons is suppressed, but while these excitons diffuse 
out of the optically-pumped area, they cool down and 
eventually recombine resulting in a local increase of the 
photoluminescence signal [^. In other words, those ex¬ 
periments did not involve Bose-Einstein condensation of 
excitons. Formation of the stationary exciton conden¬ 
sate cloud with the size comparable with the size of the 
laser pumping spot was treated in the framework of the 
Thomas-Fermi approximation in Refs. [13 and Ring 
polariton patterns were also excited by an elliptic ring 
laser spot in experiments Ref. [sj. In contrast to all 
mentioned cases, the traveling soliton-like condensates 
reported in this paper have a characteristic size smaller 
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than the excitation spot and are not caused by the ki¬ 
netics of the free carriers or by the interactions with the 
structural inhomogeneities. We infer that the soliton for¬ 
mation is caused by the interplay of the Bose-stimulated 
growth of the condensate that tends to increase the con¬ 
densate particle density and the exciton-exciton repulsion 
that leads to the particle spreading in the system. 

It is known that the condensation in infinite two- 
dimensional systems is impossible due to the phase fluc¬ 
tuations that destroy the long-range order in the system 
[^ . To achieve the exciton condensation in quasi-two- 
dimensional QW structures and to reach high enough 
exciton densities, at which the system undergoes the 
BEC transition, most of experiments were conducted 
with excitons localized in an external trapping poten¬ 
tial. A few techniques were utilized to create the exciton 
traps, including the application of non-uninform mechan¬ 
ical stress to a semiconductor sample [1^ or electrostatic 
traps [ 2 ^ [ 7 ^. In these cases, the shape of the exciton 
trap varied from a radiallv-symmetric one to a lozenge¬ 
shaped, elongated trap . 


IV. CONCLUSIONS 

To summarize, a trapped non-equilibrium condensate 
can exhibit Bose-stimulated self-organization under the 
conditions where it is driven by a homogenous external 
sources in a range of frequencies. This self-organization 
results in formation of spatially-localized structures such 
as solitary condensate waves - humps or rings. These 
solitons oscillate in a trap in a manner similar to a clas¬ 
sical particle and also can spontaneously rotate. If the 
asymmetry of the trap exceeds e > 0.5 a conventional 
fluctuating elongated dipolar exciton cloud is formed in¬ 
stead of a solitary wave. The dynamics of the conden¬ 


sate also depends on the pumping rate. It was found that 
at increased particle creation rates, the soliton-like con¬ 
densate dynamics is changed to a turbulent regime. The 
transition to turbulence at high rates can be attributed 
to the excitation of multiple interacting degrees of free¬ 
dom in the system due to nonlinearity. Turbulent exciton 
condensates have earlier been considered in Refs. [U and 

The spontaneous rotation of the soliton condensate ob¬ 
served in our simulations can be detected in experiments 
via polarization of the emitted light due to the electron- 
hole recombination photoluminescence Bm. Recently, 
polarized light has been applied to the multiplexin g in - 
formation transfer through an optical waveguide [771 . l78|| . 
The future work on the self-organized dipolar exciton 
condensate might elucidate the ways to control the soli¬ 
ton rotation, to be utilized as a source of coherent po¬ 
larized light with a given angular momentum m. It 
is of interest for prospective applications in the above- 
mentioned information technology to investigate how the 
self-organized condensate reacts on the pulsed or modu¬ 
lated pumping, to create a controllable source of coherent 
polarized light signals. Self-organized soliton condensate 
patterns can also potentially be applicable for the infor¬ 
mation transfer in exciton-based integrated circuits [t^. 
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